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Abstract. The unitary isobar model MAID2007 has been developed to analyze the world data of pion 
photo- and electroproduction. The model contains both a common background and several resonance terms. 
The background is unitarized according to the K-matrix prescription, and the 13 four-star resonances with 
masses below 2 GeV are described by appropriately unitarized Breit-Wigner forms. The data have been 
analyzed by both single-energy and global fits, and the transverse and longitudinal helicity amplitudes have 
been extracted for the four-star resonances below 2 GeV. Because of its inherent simplicity, MAID2007 is 
well adopted for predictions and analysis of the observables in pion photo- and electroproduction. 
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1 Introduction 

Our knowledge about the excitation spectrum of the 
nucleon was originally provided by elastic pion-nucleon 
scattering pQ. All the resonances listed in the Particle 
Data Tables [2] have been identified by partial-wave 
analyses of this process with both Breit-Wigner and pole 
extraction techniques. From such analyses we know the 
resonance masses, widths, and branching ratios into the 
ttN and 7T7tA^ channels. These are reliable parameters 
for the resonances in the 3- and 4-star tiers, with only 
' few exceptions. In particular, there remains some doubt 
about the structure of two prominent resonances, the 
Roper Pn(1440), which appears unusually broad, and 
the Sii(1535), where the pole can not be uniquely 
determined, because it lies close to the rjN threshold. 

On the basis of these relatively firm grounds, addi- 
tional information can be obtained for the electromagnetic 
(e.m.) <yNN* couplings through pion photo- and elec- 
troproduction. These couplings are described by electric, 
magnetic, and charge transition form factors, G* E (Q 2 ), 
G* M (Q 2 ), and Gq(Q 2 ), or by linear combinations thereof 
as helicity amplitudes Ai/ 2 (Q 2 ), A 3 / 2 (Q 2 ), and 5*1/2 (Q 2 )- 
So far we have some reasonable knowledge of the trans- 
verse amplitudes A\j 2 and A 3 / 2 at the real photon point, 
which are tabulated in the Particle Data Tables. For finite 
Q 2 the information found in the literature is scarce and 
until recently practically nonexistent for the longitudinal 
amplitudes Si/ 2 . But even for the transverse amplitudes 
only few results have remained firm over the recent 
years, such as the G* M form factor of the P33(1232) or 
Z\(1232) resonance up to Q 2 » 10 GeV 2 , the A 1/2 (Q 2 ) 
for the 5n(1535) up to Q 2 w 5 GeV 2 , and the helicity 



asymmetry A(Q 2 ) for the resonances Z?i3(1520) and 
F 15 (1680) up to Q 2 « 3 GeV 2 [3J. Frequently also data 
points for other resonances, e.g., the Roper resonance, are 
shown together with quark model calculations. However, 
the statistical errors are often quite large and the model 
dependence of the analysis may be even larger. In this 
context it is worth mentioning that also the notion of a 
'data point' is somewhat misleading because the photon 
couplings and amplitudes can only be derived indirectly 
by a partial- wave analysis. It is in fact prerequisite to 
analyze a particular experiment within a framework 
based on the "world data" . The only exception from this 
caveat is the Z\(1232) resonance. For this lowest-lying 
and strongest resonance of the nucleon, the analysis is 
facilitated by two important constraints: the validity of 
(I) the Watson theorem at the 1 % level and (II) the 
truncation of the multipole series to S and P waves as a 
good first-order approximation. With these assumptions 
the e.m. couplings have been directly determined in 
the real photon limit by a complete experiment with 
polarized photons and detecting both neutral and charged 
pions in the final state, thus allowing also for an isospin 
separation [4]. Moreover, a nearly complete separation 
of the possible polarization observables has recently 
provided the basis to extend such a "model-independent" 
analysis also to electroproduction [5]. However, we are 
still far from such a situation for all the higher reso- 
nances. Neither are the mentioned constraints valid nor 
are we close to a complete experiment. Until recently 
the data base was rather limited, the error bars were 
large, and no data were available from target or recoil 
polarization experiments. Even now there exist only very 
few data points from double-polarization experiments at 
energies above the Z\(1232). However, the situation for 
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unpolarized e + p — > e' + p + tt° reaction has considerably 
improved, mainly by new JLab experiments in all three 
halls A, B, and C. These data cover a large energy range 
from the Z\(1232) up to the third resonance region with 
a wide angular range in 9^. Furthermore, electron beam 
polarization has been used in several experiments at 
JLab, MIT/Bates, and MAMI/Mainz. Because of the 
large coverage in the azimuthal angle by the modern 
large-acceptance detectors, a separation of all 4 partial 
cross sections in the unpolarized experiment becomes 
possible. But even without a Rosenbluth separation of 
the transverse (ctt) and longitudinal (ol) cross sections, 
there is an enhanced sensitivity to the longitudinal 
amplitudes due to the interference terms cflt and olt 1 - 
Such data are the basis of our new partial-wave analysis 
with an improved version of the Mainz unitary isobar 
model MAID. 

We proceed by presenting a brief history of the uni- 
tary isobar model in Sect. 2. The formalism of pion photo- 
and electroproduction is summarized in Sect. 3. In the fol- 
lowing Sect. 4 we present our results for photoproduction 
as obtained from the latest version MAID2007, and in 
Sect. 5 this analysis is extended to electroproduction. We 
conclude with a short summary in Sect. 6. 



2 History of MAID 

MAID98 

In 1998 the first version of the Unitary Isobar 
Model was developed and implemented on the web 
to give an easy access for the community. MAID98 
was constructed with a limited set of nucleon reso- 
nances described by Breit-Wigner forms and a non- 
resonant background constructed from Born terms 
and t-channel vector-meson contributions [6]. In or- 
der to have the right threshold behavior and a rea- 
sonable description at the higher energies, the Born 
terms were introduced with an energy-dependent mix- 
ing of pseudovector and pseudoscalar ttNN coupling. 
Each partial wave was unitarizcd up to the two- 
pion threshold by use of Watson's theorem. Specif- 
ically, the unitarization was achieved by introduc- 
ing additional phases (j)R in the resonance ampli- 
tudes in order to adjust the phase of the total am- 
plitude. Only the following 4-star resonances were in- 
cluded: P 33 (1232), P n (1440), £> 13 (1520), 511(1535), 
Sii(1650), Fi 5 (1680), and L> 33 (1700). The e.m. ver- 
tices of these resonances were extracted from a best 
fit to the VPI/GWU partial- wave analysis [7]. For the 
P 33 (1232) resonance, we determined the following ra- 
tios of transition amplitudes: (I) electric quadrupole to 
magnetic dipole transition, Rem = E2/M1 = —2.2%, 
and (II) electric Coulomb to magnetic dipole transi- 
tion, R C m = C2/MI = -3.6%, independent of the 4- 
momentum transfer Q 2 . The Q 2 -dependence of the res- 
onance amplitudes in the second and third resonance 
regions was expressed in terms of the quark electric 



and magnetic multipoles [8] . The non-unitarized back- 
ground contributions were determined using standard 
Born terms and vector-meson exchange. In order to 
preserve gauge invariance, the Born terms were ex- 
pressed by the usual dipole form for the Sachs form 
factors, and both the pion and the axial form fac- 
tor were set equal to the isovector Dirac form factor, 
FAQ 2 ) = Ga(Q 2 ) = Ff(Q 2 ) - F?(Q 2 ). 
MAID2000 

In this version of MAID, the background contribution 
was unitarized for the multipoles up to F waves ac- 
cording to the prescription of K-matrix theory. The 
S'-wave multipoles Eq + and Lq + were modified in or- 
der to improve their energy dependence in the thresh- 
old region. With the new unitarization procedure, the 
pion photoproduction multipoles of SAID and some 
selected data for pion photo- and electroproduction in 
the energy range up to W = 1.6 GeV were fitted [9]. 
The ratios of the A(1232) multipoles were found to be 
Rem — —2.2% and Rcm = —6.5%, still independent 
ofQ 2 . 

MAID2003 

In accordance with results of Ref. [10 , the Q 2 de- 
pendence of the electric and Coulomb excitations 
of the Zi(1232) resonance was modified. The ratio 
Rem was found to change sign at Q 2 « 3.3 GeV 2 
from negative to positive values, and Rcm decreased 
from -6.5% at Q 2 = to -13.5% at Q 2 = 
4 GeV 2 . Moreover, the following 4-star resonances 
were included in MAID2003: 5 3 i(1620), L>i 5 (1675), 
Pi 3 (1720), P 35 (1905), P 3 i(1910), and P 37 (1950). In 
contrast to previous versions, the helicity amplitudes of 
all 13 resonances were input parameters and their Q 2 
dependence was parameterized by polynomials. With 
this new version of MAID we directly analyzed all the 
pion photo- and electroproduction data available since 
1960, and for the first time we made local (single en- 
ergy) and global (energy dependent) fits, independent 
of the GWU/SAID group. 
MAID2005 

The Q 2 dependence of the Sachs form factors in the 
Born terms was replaced by the more recent param- 
eterization of Ref. [11], and at the e.m. vertices of 
the pion-pole and seagull terms, realistic pion and ax- 
ial form factors were introduced. As a result the de- 
scription of charged pion electroproduction was much 
improved. On the basis of a large amount of new 
data from MIT/Bates, ELSA/Bonn, Grenoble, Mainz, 
and Jefferson Lab, we performed new local and global 
as well as single-Q 2 fits and obtained a better de- 
scription of the data in the energy range 1.6 GeV< 
W <2 GeV |H]. 
MAID2007 

The present version of MAID is presented in some de- 
tail in the following sections. As far as pion photo- 
production is concerned, this version is identical to 
MAID2005. The main changes are related to the Q 2 
evolution of e.m. form factors. In particular, the Q 2 - 
dependence of the Z\(1232) transition form factors are 
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remodeled to be consistent with the Siegert theorem. 
As a result the ratio Rcm decreases sharply if Q 2 
approaches zero. Furthermore, our analysis of the re- 
cent high-Q 2 data [13] led to the conclusion that Rem 
remains negative in the range of the existing experi- 
ments. 



3 Formalism for pion photo- and 
electroproduction 

Let us first define the kinematics of pion photo- and elec- 
troproduction on a nucleon, 



with the background transition potential and vi?„ the 



1 *(k)+N(p)^7T(q)+N'(p'), 



(1) 



where the variables in brackets denote the 4-momenta of 
the participating particles. In the pion- nucleon center-of 
mass (cm.) system, we define 

ft* = K, k) , if = K, q) , pi* = (E N ,-k) , (2) 



(3) 



-ml, (4) 

j 1 i En the total cm. energy and Q 2 = 
k 2 - cj 2 > the square of 4-momentum deposited by the 
photon at the nucleon vertex, also referred to as the vir- 
tually of the photon. In order to simplify the notation, 
we use the following abbreviations: 




k w = k(W, 0) 



W 2 -m 2 
2W 



contribution of the "bare" resonance excitation. The re 



suiting t-matrix can be decomposed into two terms [15j 



t 7 „(W) = t*(W)+t*(W).. 



(9) 



where 



t^(W) = v^(W)+v* w (W)g (W)U N (W), (10) 
t*(W) = v*(W)+v*(W)g (W)U N (W), (11) 



The basic equations used for MAID2007 are taken from 
the dynamical Dubna-Mainz- Taipei (DMT) model [TOlfMl 
1 1 5 j . In this approach the t-matrix for pion photo- and 
electroproduction takes the form 

t^(W) = v^(W) + v^{W) g (W) U N (W) , (7) 

with w 77r the transition potential for the reaction j* N —> 
nN, tnN the ttN scattering matrix, and go the free irN 
propagator. In a resonant channel the transition potential 
u 77r consists of two terms, 



with including the contributions from both the non- 
resonant background and the j*NR vertex renormaliza- 
tion. The decomposition in resonance and background 
contributions is not unique, however, our definition has 
the advantage that all the processes starting with the e.m. 
excitation of a bare resonance are summed up in t^ v . 



3.1 Unitarized background 

The multipole decomposition of Eq. (flTj]) yields the back- 
ground contribution to the physical amplitudes in the 
channels a = ($,£, j, I) [Sj, where £, j and / denote the 
orbital momentum, the total angular momentum, and the 
isospin of the pion-nucleon final state, and £ stands for 
the magnetic (£ = M), electric (£ = E), and Coulomb or 
"scalar" (£ = S) transitions, 



&(q,k,Q 2 ) =v*' a (q,k,Q 2 ) 



(12) 



dq 



, Q ,2 t% N (q,q';W)vB i a (q',k,Q 2 ) 
W-W„ N {q') + ie 



describing the momentum of a real photon, and 

k R = k(M R , 0) , q R = q{M R ) , (6) only, i.e., 

for the real photon and pion momenta at the resonance 
position, W = Mr. 



where W^mW) is the hadronic cm. energy in the interme- 
diate state. The pion electroproduction potential vfy a is 
constructed as in Ref. [B] and contains contributions from 
the Born terms described by an energy-dependent mix- 
ing of pseudovector (PV) and pseudoscalar (PS) irNN 
coupling as well as t-channel vector meson exchange. The 
(5) quasi-potential v^ a depends on 5 parameters: the PV-PS 
mixing parameter A m as defined in Eq. (12) of Ref. [6] and 
4 coupling constants for the vector-meson exchange. The 



on-shell parts of v^ a and t^ a depend on two variables 



771 

l! ' a (q,k,i 

7T fofc,' 



t B,a 



(W,Q 2 



v B, a 

t^ a (W,Q 2 ) 



(13) 
(14) 



v yw (W) = v*(W)+vK(W). 



(8) 



The Q 2 evolution of the s- and u-channel nucleon pole 
terms of the background is described by the form factors 
of Ref. At the e.m. vertices of the pion-pole and 
seagull terms we apply a monopole form for the pion 
form factor and a dipole form for the axial form factor, 
while the standard dipole form factor is used for the 
vector-meson exchange. 

We note that the background contribution of MAID98 
was defined by i^f(MAID98) = vfy*(W,Q 2 ) and as- 
sumed to be a real and smooth function. The unitariza- 
tion of the total amplitude was then provided by an ad- 
ditional phase <p a in the resonance contribution such that 
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the phase of the total amplitude had the phase S a of the 
respective ttN scattering state. In MAID2007, however, 
the background contributions are complex functions de- 
fined according to K-matrix theory, 



with = \/3/2 and for isospin | and 5, re- 

spectively. The energy dependence of the partial width is 
given by 



t»;<*(W > Q*) = v*><»(W,Q 2 )[l + it« N (W)], (15) ^ N (W) = Ar r R (£) + (f|±f) 



Mr 
W ' 



(20) 



where the pion-nuclcon clastic scattering amplitudes, 
t" N = [r]a exp(2iS a ) — l]/2i, are described by the phase 
shifts S a and the inelasticity parameters r\ a taken from 
the GWU/SAID analysis [16]. The assumed structure of 
the background corresponds to neglecting the principal 
value integral in the pion-rescattering term of Eq. (fT2"|) . 
Our previous studies of the P-wave multipoles in the 
(3,3) channel [T7J1[T5] showed that the "pion cloud" 
contributions of the principal value integral are effectively 
included by the dressing of the ^NN* vertex. 

Furthermore, the threshold behavior of the S waves 
was improved. The results of the dynamical ap- 
proaches [17] show that the pion cloud contributions 
are very important to obtain a good description of the 
Po+ multipole in the ir°p channel. For this purpose we 
have introduced the following phenomenological term: 



A 



B 2 q 2 



G D {Q 2 ) 



(16) 



with A and B free parameters fixed by fitting the low- 
energy 7T° photoproduction data, and Gd the standard 
nuclcon dipole form factor. The threshold correction for 
the Lq+ multipole we will consider later in Sect. 5.3. As a 
result the background contribution of MAID now depends 
on 8 parameters. We furthermore account for the cusp 
effect in the tt°p channel appearing at the ir + n threshold 
by the term [IMT9] 



pcusp 
^0+ 



77T 1 



'1- 



(17) 



where cu c — 140 MeV is the ir° cm. energy at the cusp and 
o-ttN = 0A2A/m 7T + the pion charge- exchange amplitude. 



3.2 Resonance contributions 

For the resonance contributions we follow Ref. [5J and as- 
sume Breit-Wigner forms for the resonance shape, 

t»f(w,Q*) = ftQ 2 )^#|^f^«, 
' m r -w z - iMr r to t 

(18) 

where f^N (W) is the usual Breit-Wigner factor describing 
the decay of a resonance with total width r tot (W), partial 
ttN width r^N, and spin j, 



1 few m r„ 



1/2 



with Jfl = Ztot(Ma) and Xr a damping parameter and 
(3k the single-pion branching ratio. The expression for the 
total width r tot is given in Ref. [B]. The ^yNN* vertex is 
assumed to have the following dependence on W: 



f lN {W) 



kw 



x% + 4 

X 2 A- h 2 
A R + K W 



(21) 



where n is obtained from a best fit to the real photon 
data, and with the normalization condition / 7 Ar(Afff) = 1. 
The phase 4>r{W) in Eq. (fl5|) is introduced to adjust 
the total phase such that the Fermi- Watson theorem is 
fulfilled below two-pion threshold. For the S- and P-wave 
multipoles we extend this unitarization procedure up to 
W = 1400 MeV. Because of a lack of further information, 
we assume that the phases 4>R are constant at the higher 
energies. In particular we note that the phase <j)R for the 
P 33 (1232) excitation vanishes at W = M R = 1232 MeV 
for all values of Q 2 . For this multipole we may even apply 
the Fermi- Watson theorem up to W w 1600 MeV because 
the inelasticity parameter r\ a remains close to 1. For the 
D- and P-wave resonances, the phases (f)R are assumed 
to be constant and determined from the best fit. 

Whereas MAID98 [6] included only the 7 most 
important nucleon resonances, essentially with only 
transverse e.m. couplings, our present version contains 
all 13 resonances of the 4-star tier below 2 GeV with 
transverse _ electric (*4^ = Ei±)_, transverse magnetic 
(A% = Mi±), and Coulomb (A R = Si±) couplings: 
P 33 (1232), Pn(1440), Pi 3 (1520), S n (1535), S 31 (1620), 
5n(1650), Pi 5 (1675), Pi 5 (1680), P 33 (1700), Pi 3 (1720), 
P 35 (1905), P 3 i(1910), and P 37 (1950). Because we deter- 
mine the isovector amplitudes from the proton channels, 
the number of the e.m. couplings is 34 for the proton 
and 18 for the neutron channels, that is 52 parameters 
altogether. These are taken to be constant in a single-Q 2 
analysis, e.g., in photoproduction but also at any fixed Q 2 
if sufficient data are available in the chosen energy and 
angular range. Alternatively, the couplings have also been 
parameterized as functions of Q 2 , as is discussed in Sec. 5. 

The more commonly used helicity amplitudes Ax/2, 
A3/2, and S1/2 are given by linear combinations of the 
e.m. couplings A^. These relations take the form 



A e+ 



-\(£ + 2)Et + + £M e+ ], 



4/2 = 2 VW+?){Ei+ ~M i+ ), 



(22) 
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Table 1. The reduced e.m. amplitudes A a defined by Eq. (|18p in terms of the helicity amplitudes 



iV* 


E 


M 


s 


S11/S31 

P13/-P33 

Pn/Psi 
D 13 /D 33 

D 15 /D 35 

F15/F35 

F 17 /F 37 


~A/2 

5(71^3/2 - A/2) 

-i(V3A 3/2 +A 1/2 ) 
K75A/2 - A/2) 
-i(V2A/ 2 + A/ 2 ) 

KVIA/2-A/2) 


-|(^3A/2 + A /2 ) 

A/2 

-KV2A/2 + A/2) 

-3(Vl^3/2 + A/2) 


—V2S1/2 

— r A V2 

— T 6 i/2 

_ 272 Sl /2 



for resonances with total spin j 



and 



1 , 



A/2 
A3/2 



= -[(t + i)M t --(e-i)E t 



= -^y/{e-l)(£+l)(E t - + Mi-) 



(23) 



V2 



S f 



for total spin j = £— i. The inverse relations for the partial 
waves are listed in Table [1] The helicity amplitudes are 
related to matrix elements of the e.m. current between 
the nucleon and the resonance states, e.g., as obtained in 
the framework of quark models, 



A 



1/2 



A5/2 = 



5*1/2 = - 



2na, 



, u <R, 1 -\J + \N,- 1 ->C, 



2na x 



^ <R,^\J+\N,l> C, (24) 



where J+ = — -^=(J X + ij y ) and a em = 1/137. However, 
these equations define the couplings only up to a phase 
C, which in principle can be obtained from the pionic de- 
cay of the resonance calculated within the same model. 
Because this phase is ignored in most of the literature, 
the comparison of the sign is not always meaningful, espe- 
cially in critical cases such as the Roper resonance whose 
correct sign is not obvious from the data. In contrast with 
MAID98 and MAID2000, our present version uses the he- 
licity amplitudes A/2, A/2, an d £1/2 f° r photoproduction 
as input parameters, except for the P33(1232) resonance 
which is directly described by the 3 e.m. amplitudes A a . 



4.1 Data base for pion photoproduction and fit 
procedure 



The main part of the photoproduction data was taken 
from the GWU/SAID compilation of SAID2000, which 
includes the data published between 1960 and 2000, a to- 
tal of 14700 data points. A separation of these data in 
different physical channels and observables is given in Ta- 
ble [2] In the following years the data base was extended 
by including recent results from MAMI (Mainz) [20U2 111221 
GRAAL (Grenoble) [24]|25]. LEGS (Brookhaven) [26], 
and ELSA (Bonn) [27] as listed in Table El Altogether 
4976 more data points were added. As a result our full 
data base contains 19676 points within the energy range 
140 MeV< E~ <1610 MeV. 



Table 2. Number of data points from the SAID2000 data base 
for differential cross sections (da), photon asymmetries (£), 
target asymmetries (T), and recoil asymmetries (P). 



channel 


da 


E 


T 


P 


total 


rm + 


4646 


760 


645 


205 


6256 




4936 


673 


353 


540 


6502 




1554 


206 


94 


88 


1942 



Table 3. Number of data points collected after 2000 for differ- 
ential cross sections (da), photon asymmetries (£), and helicity 



asymmetries Aa — da 1 / 2 — da- 



?3/2 • 



4 Partial-wave analysis of pion 


channel 


range _B 7 (MeV) 


data points (observable) 


Ref. 


pn u 


202-790 


1129 (da) + 357 (£) 


H2 


photoproduction data 




pn 


310-780 


174 (da) + 138 (At) 


|2T] 




pil + 


180-450 


205 (da) + 129 (Aa) 


[22] 


The unitary isobar model MAID2007 has been developed 


pTT + 


463-783 


204 (da) + 102 (Aa) 


[23] 


to analyze the world data of pion photo- and electropro- 


piT + 


800-1454 


237 (£) 


[24] 


duction. In this section we fix (I) the background parame- 


p7T° 


555-1541 


861 (da) + 469 (E) 


[25] 


ters and the helicity amplitudes for pion photoproduction 


mv~ 


285-769 


300 (da) 


[26] 


(Q 2 = 0) and (II) the dependence of the resonance con- 


pn° 


513-1575 


671 (da) 




tributions on the cm. energy W. These results are then 


total 




4976 





generalized to pion electroproduction in the next section. 
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Our strategy for the data analysis is as follows. First, 
we try to find a global (energy dependent) solution by fit- 
ting all the data in the range 140 MeV< E 7 <1610 MeV. 
This allows us to determine the phase of the multipoles, 
i.e., the ratio Imt" /Ret" above the two-pion thresh- 
old. At the lower energies this phase is constrained by the 
ttN scattering phase. In a second step we perform local 
(single energy) fits to the data, in energy bins of 10 MeV 
in the range 140 MeV< £ 7 <460 MeV and of 20 MeV 
for the higher energies, by varying the absolute values of 
the multipoles but keeping the phase as previously deter- 
mined. Similar to the prescription of the SAID group we 
minimize the modified \ 2 function 

The first term on the r.h.s. of this equation is the stan- 
dard x 2 function with Ot the calculated and 6>° xp the mea- 
sured observables, 80i the statistical errors, and A^ata the 
number of data points. In the second term, AT mu it is the 
number of the varied multipoles and Xj describes the de- 
viation from the global fit. The fitting procedure starts 
with the initial value Xj — 1 corresponding to the global 
solution, and the quantity A enforces a smooth energy 
dependence of the single-energy solution. In the limit of 
A — > oo we obtain the standard x 2 td , and for A — > the 
single-energy and the global solutions become identical. 
The optimum value for A is chosen from the condition 
1 < X 2 /Xstd < 1-05- The described two-step fitting proce- 
dure can be repeated several times by adjusting the energy 
dependence of the global solution, for example by chang- 
ing the parameters Xr and n in Eqs. (|20ti21[) in order to 
improve the agreement between the global and local solu- 
tions. 



4.2 Results for pion photoproduction 

Our results for \ 2 are summarized in Table [4] by com- 
paring the local and global solutions for different energy 
ranges and channels. We recall that the number of varied 
multipoles 2V mu jt in the proton and neutron channels is 
different. Since the number of data points in the proton 
channels (jtt° and r yir + ) is about one order of magnitude 
larger than for the neutron channel (■y%~), we proceed as 
follows. First, we analyze the proton channel and extract 
the multipoles p Ej± 2 , V M\± ', Ef± 2 , and Mf± 2 as defined 
in Ref. [B]. Second, with the thus obtained values for the 

1/2 

isospin 3/2 multipoles, we extract the multipoles n E t j_ 

1/2 

and n Mi± from the neutron channel. In this way we 
minimize the pressure from the large number of proton 
data on the results in the neutron channel. The number 
of varied multipoles also depends on the energy. For 
Ery <450 MeV we vary all the S- and P-wave multipoles 

1 /2 3 /2 

plus p>n E 2 _ and E 2 _ . At the higher energies we include 
all the multipoles up to the F waves. 



Table 4. Results for \ 2 from single-energy (se) and global (gl) 
solutions. 





proton 


E 1 [MeV] 


-A mu lt 


Adata 


Xsc 


■xk 


140 — 200 


10 


990 


787 


2346 


200 — 450 


10 


5622 


5454 


14236 


4oU — ooU 


24 


6403 


6996 


22700 


850—1210 


24 


2965 


4133 


24990 


1210—1610 


24 


1454 


4737 


12594 


total 




17434 


22107 


76866 




neutron 


£ 7 [MeV] 


iVmuit 


Adata 


Xsc 


Xii 


140—200 


5 


51 


113 


151 


200—450 


5 


872 


1613 


2203 


450—850 


12 


902 


1748 


3414 


850—1210 


12 


334 


651 


2311 


1210—1610 


12 


83 


107 


583 


total 




2242 


4222 


9262 



The best fit for the background and the resonance 
parameters yields the results listed in Tables and 
IH1 respectively. The PS-PV mixing parameter and the 
vector-meson coupling constants are defined as in Ref. [6] . 
However we note that in the present version we do 
not use form factors at the hadronic vertices involving 
vector-meson exchange. In Tables [7] and [8] we corn- 
Table 5. Masses and coupling constants for vector mesons, 
PS-PV mixing parameter A m , and parameter A for the low- 
energy correction of Eq. (|16l) . 





m v [MeV] 


\v 


gvi 


gv2/gvi 




783 


0.314 


16.3 


-0.94 


P 


770 


0.103 


1.8 


12.7 


A m 


= 423 MeV 


A = 1.9 x 10" 




B = 0.71/m 



pare the helicity amplitudes obtained from MAID2003 
and MAID2007 with the results of the PDG and 
GWU/SAID [281129] analysis. As is very typical for a 
global analysis with about 20,000 data points fitted to 
a small set of 20-30 parameters, the fit errors appear 
unrealistically small. However, one should realize that 
these errors only reflect the statistical uncertainty of the 
experimental error, whereas the model uncertainty can 
be larger by an order of magnitude. We therefore do not 
list our fit errors, which in fact are very similar in the 
GW02 or GW06 fits of the SAID group [25P5]. The only 
realistic error estimate is obtained by comparing different 
analysis, such as SAID, MAID, and coupled-channels 
approaches. 

Next we present our results for the multipoles starting 
with the threshold region. In Fig. [1] we demonstrate 
the effects of the low-energy correction and the cusp 
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Table 6. Resonance masses Mr, widths Fr, single-pion branching ratios j3^, and angles 4>R as wel l as the parameters Xr,ue, 
and riM of the vertex function Eq. pill . 





Af fl [MeV] 


Pfl[MeV] 


P„ 


<Mdeg] 


X R [MeV] 


proton 

he 


n M 


neutron 

tie 


n M 


P 33 (1232) 


1232 


130 


1.0 


0.0 


570 


-1 


2 


-1 


2 


Pn(1440) 


1440 


350 


0.70 


-15 


470 









-1 


£>is(1520) 


1530 


130 


0.60 


32 


500 


3 


4 


7 


2 


Six (1535) 


1535 


100 


0.40 


8.2 


500 


2 




2 




531(1620) 


1620 


150 


0.25 


23 


470 


5 




5 




5n(1650) 


1690 


100 


0.85 


7.0 


500 


4 




4 




£>i 8 (1675) 


1675 


150 


0.45 


20 


500 


3 


5 


3 


4 


Fi 5 (1680) 


1680 


135 


0.70 


10 


500 


3 


3 


2 


2 


D 33 (1700) 


1740 


450 


0.15 


61 


700 


4 


5 


4 


5 


Pis (1720) 


1740 


250 


0.20 


0.0 


500 


3 


3 


3 


3 


P 3 5(1905) 


1905 


350 


0.10 


40 


500 


4 


5 


4 


5 


Psi(1910) 


1910 


250 


0.25 


35 


500 




1 




1 


F 37 (1950) 


1945 


280 


0.40 


30 


500 


6 


6 


6 


6 



Table 7. Proton helicity amplitudes at Q 2 — for the major 
nucleon resonances, in units 10" 3 GeV" 1/2 . The results with 
MAID2003 and MAID2007 are compared to the PDG [2 and 
GWU/SAID [29] analysis. 





PDG 


GW06 


2003 


2007 


P 33 (1232) 


A 1/2 


-135±6 


-139.1±3.6 


-140 


-140 




A3/2 


-250±8 


-257.6±4.6 


-265 


-265 


E2/M1 (%) 


-2.5±0.5 




-2.2 


-2.2 


Pii(1440) 


A 1/2 


-65 ±4 


-50.6 ±1.9 


-77 


-61 


Di 3 (1520) 


A 1/2 


-24 ±9 


-28.0 ±1.9 


-30 


-27 




A3/2 


166 ±5 


143.1 ±2.0 


166 


161 


5 U (1535) 


A 1/2 


90 ±30 


91.0 ±2.2 


73 


66 


#51 (1620) 


A 1/2 


27 ±11 


49.6±2.2 


71 


66 


Six (1650) 


A 1/2 


53±16 


22.2 ±7.2 


32 


33 


Di 5 (1675) 


A 1/2 


19 ±8 


18.0 ±2.3 


23 


15 




A3/2 


15 ±9 


21.2 ±1.4 


24 


22 


Fib (1680) 


A 1/2 


-15 ±6 


-17.3 ±1.4 


-25 


-25 




A3/ 2 


133 ±12 


133.6 ±1.6 


134 


134 


£> 33 (1700) 


A 1/2 


104 ±15 


125.4 ±3.0 


135 


226 




Aa/ 2 


85 ±22 


105.0 ±3.2 


213 


210 


Pis (1720) 


A 1/2 


18 ±30 


96.6 ±3.4 


55 


73 




Az/ 2 


-19 ±20 


-39.0 ±3.2 


-32 


-11 


F 35 (1905) 


A 1/2 


26 ±11 


21.3 ±3.6 


14 


18 




Az/ 2 


-45 ±20 


-45.6±4.7 


-22 


-28 


F 3 7(1950) 


A 1/2 


-76 ±12 




-78 


-94 




As/ 2 


-97 ±10 




-101 


-121 



effect for tt° photoproduction, as described by Eqs. (|16p 
and (fTT). respectively. The prediction of MAID98 for 
7T° photoproduction at threshold (dotted lines) lies sub- 
stantially below the data. In accordance with Ref. [T7] . 
the phenomenological term £^ rr simulates the pion 
off-shell rescattering or pion-loop contributions of ChPT. 
The cusp term of Eq. (fT7|) describes the strong energy 
dependence near n + threshold, which has its origin in 
the pion mass difference and the strong coupling with 



Table 8. Neutron helicity amplitudes at Q 2 = for the major 
nucleon resonances. GW02 are the results GWU/SAID analy- 
sis [28] . Further notation as in Tab. [7] 





PDG 


GW02 


2003 


2007 


Pn(1440) 


A 1/2 


40±10 


47±5 


52 


54 


Di 3 (1520) 




-59 ±9 


-67 ±4 


-85 


-77 




^3/2 


-139 ±11 


-112 ±3 


-148 


-154 


511 (1535) 


A 1/2 


-46±27 


-16±5 


-42 


-51 


Sxx(1650) 


A 1/2 


-15±21 


-28±4 


27 


9 


£>is(1675) 


A 1/2 


-43 ±12 


-50 ±4 


-61 


-62 




A-A/ 2 


-58 ±13 


-71 ±5 


-74 


-84 


Fi 5 (1680) 


A 1/2 


29 ±10 


29 ±6 


25 


28 




A 3/2 


-33 ±9 


-58 ±9 


-35 


-38 


Pis (1720) 


A 1/2 


1 ±15 




17 


-3 




A- i/2 


-29 ±61 




-75 


-31 



the ir + n channel. The figure shows that the off-shell 
pion rescattering substantially improves the agreement 
with the data. However, one problem still remains in the 
threshold region. The experimental photon asymmetry S 
in 7T° photoproduction at w 160 MeV takes positive 
values, whereas the MAID results are negative in this 
region. As has been demonstrated in Refs. [T7ll32j. 
this observable is very sensitive to the Mi_ multipole 
which strongly depends on the details of the low-energy 
behavior of Roper resonance, vector meson and off-shell 
pion rescattering contributions. Therefore, a slight modi- 
fication of one or all of these mechanisms can drastically 
change the photon asymmetry. 

Figures [2j4] display the results for the most important 
S and P waves in the Z\(1232) region. However, a look 

i - 1 1 /2 

at Fig. [5] shows that also the D-wave amplitudes P E 2 '_ , 

3/2 1/2. 

E' 2 _ , and n E 2 '_ , give sizable contributions in this region, 
in particular through their real parts. In these figures, we 
present the MAID and SAID global (energy dependent) 
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solutions, together with our local (single energy) fit 
obtained for energy bins of 10 MeV. In general the MAID 
and SAID results are close, which is not too surprising 
because the phases are constrained by the Fermi- Watson 
theorem. However, there are much larger differences in 

1/2 1/2 

the p E 0+ and p E 2 _ amplitudes, which indicates that the 
present data base is still too limited to determine these 
background amplitudes in a reliable way. 

More substantial discrepancies between the MAID and 
SAID analyses are found in the second and third resonance 
regions. A detailed comparison of the two models is shown 
in Figs.[6l fT3l As pointed out in Sect. 3.1, it is prerequisite 
to know the phases of the multipoles in order to get correct 
single-energy solutions above the two-pion threshold. In 
MAID2007 these phases are determined by Eqs. ©, (|T5|l . 
and (fl8]) . The SAID analysis is based on the following 
parametrization of the partial wave amplitudes: 



tipole at E~. 



1070 MeV or W 



1700 MeV. How- 



(Born + A) (1 + it vN ) + B t vN 
+ (C + iD) {Imt nN - | Un 12 1 



(26) 



where A, B C and D are polynomials in the energy 
with real coefficients, and t^N is the pion-nucleon elastic 
scattering amplitude of Eq. (fT5"|) . As seen in the first 
resonance region, the most serious differences between 
MAID and SAID are again found for the real parts of the 
multipoles P E^ and P E\^ . We have checked the phases 
of these multipoles by independent calculations on the 
basis of dispersion relations 33 , 34J . The result confirmed 
our phase relations. Concerning the small amplitudes, the 
most sizable differences between SAID and MAID are in 

3/2 1 /2 1/2 

the M x '_ pMjl and P E^ multipoles. In the neutron 
channel, the largest differences are in the multipoles 

rt-Eo+ 2 , nE\J\ nE\+ , and n M\^ . In the last two cases 
this is due to the large contribution from the Pi3(1720) 
resonance which is not found in the SAID analysis (see 
Table [3). 

Let us finally discuss the possible contributions of the 
weaker resonances. As discussed in Ref. [35], the two addi- 
tional S\\ resonances found with masses of about 1800 and 
2000 MeV might also show up in pion photoproduction. 
This conclusion was mainly based on the single-energy so- 
lution of the SAID group. As illustrated by Fig. QjJ our 
present analysis requires only one additional S\\ resonance 
with mass Mr 1950 MeV, and our single-energy solu- 
tion shows no resonance at Mr w 1800 MeV. Of course, 
the solution of the problem is certainly correlated with 
the way how the resonance and background contributions 
are separated. As demonstrated both in Ref. [35] and by 
Fig. [6] the background is very important for this particu- 
lar channel. In this context we recall that we use the same 
form of the unitarized background contribution (Born, u> 
and p exchange) for all the partial waves. Another in- 
teresting topic deserving further experimental and the- 
oretical studies, concerns the Roper or Pu channel. As 
clearly seen in Fig. [151 both our and the SAID analy- 

1 /2 

sis yield a second resonance structure of the p M 1 _ mul- 



ever, our analysis yields a very small width of /tot ~ 30- 
60 MeV, whereas the PDG lists r tot ~ 50-250 MeV for 
this resonance. Moreover, also the helicity amplitude dif- 



fers, whereas our result is A 1 



I 2 



-0.024 GeV" 1 / 2 , the 
PDG lists 0.009 ± 0.022 GeV" 1 / 2 . Of course, these num- 
bers do strongly depend on the values for the single-pion 
branching ratio. On the other hand, we do not anticipate 
large effects from different definitions of the background 
in this channel, because the background contribution is 
very small in the resonance region (see Fig. [7]). 



5 Partial-wave analysis of pion 
electroproduction 

In most of the pion electroproduction experiments the five- 
fold differential cross section was measured. However, dif- 
ferent conventions exist for the partial cross sections, and 
therefore we recall the definitions used in MAID. For an 
unpolarized target the cross sections written as the prod- 
uct of the virtual-photon flux factor E v and the virtual 
photon cross section da v /dfi n [36] . 



da 



da,, 



dn) dEj dQ« df2 



da v 
df}.,. 



dan 



da i 



dH-,, dQ 
daTT 



+ y/2e(l + e) 



da 



LT 



dfl^ 



(27) 



cos^ (28) 



dQ^, 



cos 23V + hy/2e(l-e) ■ sin 4> n . 



where e and h describe the polarizations of the virtual 
photon and the electron, respectively. We further note 
that the hadronic kinematics is expressed in the cm. sys- 
tem, whereas the electron and virtual photon kinematics 
is written in the lab frame, as indicated by L in the fol- 
lowing variables: the initial and final electron energies Ef 
and Ej, respectively, the electron scattering angle 6l, the 
photon energy u>l — E\ — Ej, and the photon three- 
momentum k^. With these definitions the virtual photon 
flux and the transverse photon polarization take the form 



1 



i 



i 



J Q- 



■ tan z ^ 



2tt 2 El Q 2 1 - e 



(29) 



As in our previous notation [36], the flux is denoted 
by the photon "equivalent energy" in the lab frame, 
K = Kh = (W 2 — m 2 )/2m as originally introduced by 
Hand |37j . Another definition was given by Gilman [38 
who used K = Kq = | | . 



The first two terms on the r.h.s. of Eq. ([2"5| are the 
transverse (T) and longitudinal (L) cross sections. They 
do not depend on the pion azimuthal angle ^ . The third 
and fifth terms describe longitudinal-transverse interfer- 
ences (LT 7 LT'). They contain an explicit factor sin^ and 
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therefore are vanishing along the axis of momentum trans- 
fer. The same is true for the fourth term, a transverse- 
transverse interference (TT) proportional to sin 2 9^. It is 
useful to express these 5 cross sections in terms of hadronic 
response functions depending only on 3 independent vari- 
ables, i.e., Ri — Ri(Q 2 ,W,0 n ). The corresponding rela- 
tions take the form 

d(JT _ _q_ R ^ d&TT _ _q_ ^ da L _ q Q 2 R ^ 
dfl-x kw ' dOn kw ' dfi^ kw ^ 2 ' 

dcJLT q Q D da LT ' q Q , . 

~7cT = I Klt ' ac = I Klt ' ■ ^0) 

As a result of this equation, the longitudinal (L) and 
longitudinal-transverse (LT and LT 1 ) response functions 
must be proportional to w 2 and w 7 , respectively, in order 
to avoid non-physical singularities at the energy for which 
the cm. virtual photon energy passes through zero. The 
5 response functions may be expressed in terms of 6 
independent CGLN amplitudes F%, Fq [39], or in terms 
of the helicity amplitudes Hi , . . . , Hq , which are linear 
combinations of the CGLN amplitudes. The relevant 
expressions can be found in Refs. [BlflPIO] . 



Table 9. The number of data points, Adata, and the \ 2 value 
per data point obtained with MAID2003 and MAID2007. 



Ref. 
channel 


W (MeV) 
Q 2 (GeV 2 ) 


Adata 

observables 


XV^data (2003) 
xV^data (2007) 


SAIDOO 


1074-1895 


13152 


3.238 


P7T 


0.1-4.3 


da, ... 


3.172 


SAIDOO 


1125-1975 


5464 


3.297 


nir + 


0.117-4.4 


da, ... 


4.188 


Bonn02 [42] 


1153-1312 


4914 


1.378 


P7T 


0.63 


da 


1.400 


CLAS02 03] 


1110-1680 


31810 


1.907 


pn° 


0.4-1.80 


da 


1.952 


CLAS03 Si] 


1100-1660 


223 


4.881 


pn° 


0.4-0.65 


da LT , 


3.490 


CLAS04 [45] 


1100-1660 


224 


4.879 


mv + 


0.4-0.65 


da LT i 


2.196 


CLAS06 g6] 


1110-1570 


4179 


10.04 


nir + 


0.3-0.60 


da 


4.954 


CLAS06 [13] 


1110-1390 


8491 


1.691 


pn° 


3.0-6.0 


da 


1.335 


total 


1074-1975 


68457 


2.724 


pn°, nir + 


0.1-6.0 


da, ... 


2.437 


SAIDOO 


1253-1976 


799 


2.100 




0.54-1.36 


da 


2.264 



5.1 Data base for pion electroproduction and fit 
procedure 

The main part of our data base for pion electroproduction 
includes the compilation of the GWU/SAID group [H] in 
2000 and recent data from Bonn and JLab (see Table [9|). 
Altogether this base contains about 70000 data points 
within the energy range 1.074 GeV< W < 2 GeV and 
photon virtuality range 0.1 GeV 2 < Q 2 <6 GeV 2 . In addi- 
tion we have analyzed high precision data from Bates [471 
148] . Mainz [491150] . and JLab [5ll5T]. Our fitting proce- 
dure was as follows. In a first step we fitted the data sets 
at constant values of Q 2 (single-Q 2 fit). This procedure 
is similar to the partial-wave analysis for pion photopro- 
duction except for the additional longitudinal couplings of 
the resonances. Second, we introduced a smooth Q 2 evolu- 
tion of the e.m. transition form factors and parameterized 
the 3 helicity amplitudes accordingly. In a combined fit 
with the complete electroproduction data base and infor- 
mation from the single-Q 2 fits we finally constructed the 
Q 2 -dependent solution (super-global fit). This new solu- 
tion (MAID2007) was then compared with the previous 
solution (MAID2003) in terms of x 2 as presented in Ta- 
ble [5] In most cases the new fit improves the description 
of the data, in particular for the nn + channel. 



5.2 Results for the Z\(1232) form factors 

In the literature the e.m. properties of the NA(1232) tran- 
sition are described by either the magnetic {G* M ), electric 
(G* E ), and Coulomb {G* c ) form factors or the helicity am- 
plitudes A1/2, A3/2, and S1/2, which can be derived from 



the reduced e.m. amplitudes A a as defined by Eq. (fl8|) . It 
is worthwhile pointing out that these amplitudes are re- 
lated to the multipoles over the full energy region, that is, 
they are the primary target of the fitting procedure. The 
form factors and helicity amplitudes are then obtained by 
evaluating the reduced e.m. amplitudes at the resonance 
position W — A/,4^1232 MeV. The respective relations 
take the following form: 

G* M (Q 2 ) = -ca(A 1/2 + V3A 3/2 ) = 2c A Am (Ma, Q 2 ) , 
G* E (Q 2 ) = c A (A 1/2 - 7^3/2) = -2c A A%(M A ,Q 2 ) , 

G* C (Q 2 ) = V2c A ^S 1/2 = -2c^ ^Af(M A ,Q 2 ) , 
k A k A 



with c A 



m 3 k& 



4wa em M A k 2 A 



1/2 



(31) 



and where k A = k A (Q 2 ) = k(M A ,Q 2 ) and k w = 
k(M A , 0) are the virtual photon momentum and the pho- 
ton equivalent energy at resonance. Because the Z\(1232) 
is very close to an ideal resonance, the real part of the 
amplitudes vanishes for W = M A and the form factors 
can be directly expressed by the imaginary parts of the 
corresponding multipoles at the resonance position, 



G* M (Q 2 ) - b A Im{M^ 2) (M A , Q 2 )} , 
G* E (Q 2 ) = -b A lm{E[ 3 _l 2 \M A ,Q 2 )}, 

Gh(Q 2 ) = -b A ^lm{ S [f\M A ,Q 2 )}, 



(32) 
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where 



8 to 2 £a 
3 a 



1/2 



and with — 115 MeV and qA = q(Ma) the pion mo- 
mentum at resonance. The above definition of the form 
factors is due to Ash (52]. The form factors of Jones and 
Scadron }53] are obtained by multiplying our form fac- 
tors with y / l + Q 2 /(M N + M A ) 2 - We note that the form 
factor G* c differs from our previous work [SI] by the fac- 
tor 2M A /kA in Eq. (J32J) ■ With these definitions all 3 
transition form factors remain finite at pseudo-threshold 
(Siegert limit). In the literature, the following ratios of 
multipoles have been defined: 



Rem — — 
Rsm 



CA 



M/2 - 775^3/2 



kA G* c _ V%Si/2 



2M A G* m A 1/2 + V3A 



(33) 
(34) 



i3/2 



In MAID2003 the Q 2 dependence of the e.m. ampli- 
tudes was parameterized as follows: 

At(W,Q 2 ) = A° a (l + a Q 2n ^e-^ 2 G D (Q 2 ), (35) 

where G D (Q 2 ) = 1/(1 + Q 2 /0.71GeV 2 ) 2 is the dipole 
form factor. MAID2007 follows this parametrization for 
the magnetic and electric amplitudes, although with some- 
what different values of the parameters (see Table [TO]) . 
In order to fulfill the Siegert theorem, we have however 
changed the description of the Coulomb amplitude as 
specified below. The results of MAID2003 and MAID2007 



Table 10. Parameters for the NA amplitudes given by 
Eqs. ([35]) and (02J. The amplitudes A° a are in units 10 -3 
GeV -1 / 2 , the parameters /3 and 7 in GeV -2 . For the Coulomb 
amplitude in MAID2007 we use Eq. J42J with d=4.9. 





ME S 


model 




300 -6.50 -19.50 
300 -6.37 -12.40 


2003 
2007 




-0.306 0.017 
0.01 -0.021 0.12 


2003 
2007 


7c 


0.21 0.21 0.21 
0.23 0.16 0.23 


2003 
2007 


n a 


1 1 3 
1 1 


2003 
2007 



-Tzi=130 MeV, which is different from the usually assumed 
value of about 115 MeV. Therefore, in order to compare 
with form factor values of other analyses, we scale our pre- 
dicted form factor with -y/115/130. As shown in Fig. [TH] 
G* M (Q)/i will then take the usual value of 1 to an accuracy 
of 1%. From this number we can determine the N — > A 
magnetic transition moment, fiNA — 3.46 ± 0.03, in units 
of the nuclear magneton. 



5.3 Siegert theorem and ratios Rem and R S m 

Let us next discuss our results for the Rem and Rsm 
ratios. In all previous solutions these ratios were nearly 
constant for Q 2 < 1 GeV 2 . However, calculations in ef- 
fective field theories 55,56 and dynamical models 10, 
IT51I57] indicated a rapid rise of Rsm for Q 2 — > +0. This 
dependence is rather model- independent, because it re- 
flects the behavior of the multipoles at physical threshold 
(pion momentum q — * 0) and pseudothreshold (Siegert 
limit, photon momentum k — > 0) [55] . The longitudinal 
and Coulomb multipoles are related by gauge invariance, 
k • J = LU 7 p, which leads to 



\MLl ± (W,Q 2 )=Lu 7 S I e± (W,Q 2 ) 



(36) 



Since the photon cm. energy w 7 vanishes for Q 2 — Qq = 
W 2 — m 2 , the longitudinal multipolc must have a zero 
at that momentum transfer, Lj,(W,Qo) = 0. Further- 
more, gauge invariance implies that the longitudinal and 
Coulomb multipoles take the same value in the real pho- 
ton limit, L{ ± {W,Q 2 = 0) = Sl ± {W,Q 2 = 0). Finally, 
the multipoles obey the following model-independent re- 
lations at physical threshold (q — > 0) and pseudothreshold 
(k-0): 



(eI + ,lL) — ► k l q l (l>0) 



(M/ + ,A//_)-fcV 
{El,Ll_)^kl- 2 



(*>1) 
(£=1) 

f {I > 2) 



(37) 



According to Eq. (|36[) the Coulomb amplitudes acquire an 
additional factor k at pseudothreshold, i.e., ~ kL^. 
This limit is reached at Q 2 = Q 2 t = —(W — to) 2 (pseudo- 
threshold), and because no direction is defined for k = 
0, the electric and longitudinal multipoles are no longer 
independent at this point, 



are compared in Fig. M Because our single- In the case of the NA mu ltipoles, Eq 



for G* M (Q 2 ) 

Q 2 analysis follows the global fit closely, it is not shown 
in the figure. We find an excellent agreement with the 
data, which also include the new high-Q 2 data of the 
JLab/CLAS Collaboration [13]. At this point a word of 
caution is in order. Because the form factors are extracted 
from the multipoles by Eq. ([3"2"]l . they are proportional 
to \fTA~. The MAID fit to the experimental data yields 



E t+/ Ll t+ -> 1 an d El-lH- -> ~t/(t - 1) if fc — »> 0. 

(38) 

yields the fol- 



lowing relation in the limit k — > 0: L\^ 2 ~ > E±+ —> 0{k) 

and consequently S\+ = kE^/ujj — > 0(k 2 ). Although 
the pseudo-threshold is reached at the unphysical point 
Q 2 t = -(M A - m) 2 fa -0.084 GeV 2 , it still influences 
the multipoles near Q 2 = because of the relatively small 
excitation energy of the Z\(1232). In particular we get the 



3/2 
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following relation for Q 2 



Q 2 



Rs 



M 



(3/2) 



k E- 



(3/2) 
1+ 



M 



(3/2) 
1+ 



h M{ 



(3/2) 



M A 



-R 



EM ■ 



(39) 



With increasing value of Q 2 , the Siegert relation fails to 
describe the experimental data. Moreover, it contains a 
singularity at w 7 =0, which occurs in Z\(1232) elcctropro- 
duction already at Q 2 = 0.64 GeV 2 . However, we obtain 
a good overall description by using the idea of Ref. [55] 
that the ratio Rsm is related to the (elastic) form factors 
of the neutron, 



Rsm{Q 2 ) 



mk A (Q 2 )G n E (Q 2 
2Q 2 Gl I (Q 2 ) 



(40) 



This relation gives the necessary proportionality to the 
photon momentum at small Q 2 , describes the experimen- 
tal value of the ratio over a wide range of Q 2 , and yields an 
asymptotic behavior consistent with the prediction of per- 
turbative QCD that Rsm should approach a constant for 
Q 2 — > co. This leads to the following simple parametriza- 
tion: 

with t = Q 2 1 (Am 2 ), and the parameters a and d to be 
determined by a fit to the data. On the basis of this ansatz, 
the Coulomb coupling has been modified as follows: 



Af(W,<?) = A» s : 



1 + dr kw k 



G D (Q 2 ),(42) 

w 



with parameters given in Table [T01 This leads to the mul- 
tipolc ratio 

A% 1 (1 + I3 S Q 2 \ k A 



Rsm(Q') = -rf —j- LI" hj. (43) 
A U M l + dr \l + I3 M Q Z J k? 



nv 



By construction this ratio vanishes in the Siegert limit, 
Q 2 ~ * QptJ an d approaches a (negative) constant for 
Q 2 — ► co in agreement with perturbative QCD. However, 
a word of caution has to be added at this point. The 
polynomials and gaussians used to fit the data in the 
range of low and intermediate virtualities, Q 2 < 10 GeV 2 , 
should not be expected to yield realistic extrapolations 
to the higher values of Q 2 . 

The correct Siegert limit is even more important for 
pion S-w&ve production in the threshold region, in which 
case the pseudo-threshold comes as close as Q 2 t = —to 2 w 
—0.02 GeV 2 . The term describing the pion cloud contri- 
bution has therefore been parameterized as follows: 

L%?(W, Q 2 ) = ^ e-W 2 -^*) E^(W, Q 2 ) , (44) 

where w 2 t — -Q 2 t — (W — to) 2 . From a fit to 7r° 
electroproduction data near threshold [ST], we obtain 
P = 10 GeV -2 . In the future we intend to study pion 



electroproduction near threshold in more detail. 

Figures [Tfl and PT8l display the super-global solutions of 
MAID2003 (dashed lines) and MAID2007 (solid lines) for 
the ratios Rem and Rsm in comparison with other analy- 
ses. Different from our previous solution, the ratio Rem of 
MAID2007 stays always below the zero line, in agreement 
with the original analysis of the data [T3"Il6"2"] and also with 
the dynamical model of Sato and Lee [57] who concluded 
that Rem remains negative and tends towards more nega- 
tive values with increasing Q 2 . This indicates that the pre- 
dicted helicity conservation at the quark level is irrelevant 
for the present experiments. We also analyzed the new 
data of Ref. 13J in the range of 3 GeV 2 < Q 2 < 6 GeV 2 
and found slightly decreasing values of Rem from our 
single-Q 2 analysis. In this analysis we varied both the A 
and the Roper multipoles. For the ratio Rsm both the 
super-global and the single-Q 2 solutions yield ratios that 
asymptotically tend to a negative constant. This result is 
in good agreement with the prediction of Ref. [59] (dash- 
dotted curve in Fig. [THJ) but disagrees with our previous 
solution and with the analysis of Ref. [T3]. As discussed 
before, the new solution has a large slope at small Q 2 as a 
consequence of the Siegert theorem. The following Fig. [19] 
displays the Q 2 dependence of the helicity amplitudes for 
the JV .4(1232) transition. Our single-Q 2 fit is in excellent 
agreement with the super-global solution, except for the 
values of S 1/2 at Q 2 = 0.4 and 0.525 GeV 2 . 



5.4 Results for the higher resonances 

Above the two-pion threshold we can no longer apply 
the two-channel unitarity and consequently the Watson 
theorem does not hold. Therefore, the background 
amplitude of the partial waves does not vanish at reso- 
nance as was the case for the Z\(1232) resonance. As an 
immediate consequence the resonance-background sep- 
aration becomes more mo del- dependent. In MAID2007 
we choose to separate the background and resonance 
contributions according to the K-matrix approximation. 
Furthermore, we recall that the absolute values of the 
helicity amplitudes are correlated with the values used 
for the total resonance width _Tr and the single-pion 
branching ratio On the experimental side, the data 
at the higher energies are no longer as abundant as 
in the A region. However, the large data set recently 
obtained by the CLAS collaboration (see Table \§§ 
enabled us to determine the transverse and longitudinal 
helicity couplings as functions of Q 2 for all the 4-star 
resonances below 1700 MeV. These data are available in 
the kinematical region of 1100 MeV < W < 1680 MeV 
and 0.4 GeV 2 < Q 2 < 1.8 GeV 2 . 

The helicity amplitudes for the Roper resonance 
Pn(1440) are shown in Fig. [201 Our latest super-global 
solution (solid lines) is in reasonable agreement with the 
single- Q 2 analysis. The figure shows a zero crossing of 
the transverse helicity amplitude at Q 2 » 0.7 GeV 2 and 
a maximum at the relatively large momentum transfer 
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Q 2 w 2.5 GeV 2 . The longitudinal Roper excitation rises 
to large values around Q 2 ps 0.5 GeV 2 and in fact pro- 
duces the strongest longitudinal amplitude we can find in 
our analysis. This answers the question raised by Li and 
Burkert [S3] whether the Roper resonance is a radially 
excited 3-quark state or a quark-gluon hybrid, because 
in the latter case the longitudinal coupling should van- 
ish completely. From the global fit we find the following 
parametrization for the Q 2 dependence of the Roper am- 
plitudes for the proton and neutron channels: 

A p 1/2 (Q 2 ) = Al'Ul - 1.22 Q 2 - 0.55 Q 8 ) e~ 151 « 2 , 



Sf /2 (Q 2 ) 



°x% (1 

1 0,n , 



40Q 2 + 1.5Q 8 )e- 175Q2 , (45) 



A1 /2 (Q Z ) =A" X 'P 2 (1 + 0.95 Q 
S{) 2 (Q 2 ) = S° 1/ n 2 (l + 2.98Q 2 



2\ g— 1.77Q 



-1.55Q 2 



(46) 

where Q 2 should be inserted in units of GeV 2 . The 
numerical values of the hclicity amplitudes for real 
photons are given in Table [IT] At Q 2 =0 the fit yields a 
large neutron value for the Coulomb amplitude S^ 2 , but 
with increasing Q 2 the proton and neutron amplitudes 
become comparable. 



JLab- Yerevan analysis, which was performed with both 
an isobar model and dispersion analysis. The following 



Table 12. The proton parameters for the higher resonances: 
a and P as denned by Eq. (|47p. in units GeV -2 , and S^ 2 j the 
longitudinal amplitude at Q 2 = 0, in units 10 -3 GeV _1//2 . The 
values for the transverse amplitudes A°i 2 3/2 are determined 
by the real photon physics and listed in Table [7] 





A 


1/2 


A 


3/2 


S1/2 


D l/2 


proton 


a 


P 


a 


P 


a 


P 




£>is(1520) 


7.77 


1.09 


0.69 


2.10 


4.19 


3.40 


-63.6 


Six (1535) 


1.61 


0.70 






23.9 


0.81 


-2.0 


S 3 i(1620) 


1.86 


2.50 






2.83 


2.00 


16.2 


Su(1650) 


1.45 


0.62 






2.88 


0.76 


-3.5 


■Die (1675) 


0.10 


2.00 


0.10 


2.00 


0.00 


0.00 


0.00 


Pis (1680) 


3.98 


1.20 


1.00 


2.22 


3.14 


1.68 


-44.0 


P 33 (1700) 


1.91 


1.77 


1.97 


2.20 


0.00 


0.00 


0.00 


Pis(1720) 


1.89 


1.55 


16.0 


1.55 


2.46 


1.55 


-53.0 



Table 11. Helicity amplitudes for the Pn(1440) resonance at 
Q 2 =0 in units 10~ 3 GeV~ 1/2 . 



Table 13. The neutron parameters for the higher resonances. 
The values for the transverse amplitudes A\j 2 3//2 are given in 
Table[8] Further notation as in Table [U 





^1/2 
proton neutron 


°l/2 

proton neutron 


Pn(1440) 


-61.4 54.1 


4.2 -41.5 



For all the higher resonances the transverse and lon- 
gitudinal helicity amplitudes are simply parameterized by 
the form 



A X (Q 2 ) = A° X (l + aQ 2 )e 



-0Q 2 



(47) 



The values of the fit parameters A° X: a and (3 are listed 
in Tables [12] and [13] In the following Fig. [21] we present 
the results for the fin (1535). As is also known from 77 
electroproduction, the transverse form factor falls off very 
slowly. At a virtuality of Q 2 ps 3 GeV 2 this resonance 
is much stronger than the Z\(1232) or the Z?i3(1520) 
and only comparable to the Roper. However, due to its 
much smaller width as compared to the Roper, the Sn 
dominates over the Roper at large Q 2 . This result is 
in agreement with the inclusive electroproduction cross 
section on the proton, which clearly shows the dominance 
of the Z\(1232) at small momentum transfer whereas at 
the larger momentum transfers the second resonance 
region takes over. 

In Fig. [22] we compare our results to those of Az- 
nauryan et al. [64] who used a similar set of the CLAS 
data in the second resonance region. Our super-global so- 
lutions (solid lines) agree generally quite well with the 



neutron 


A 1/2 
a P 


^3/2 
a /3 


Sl/2 

a p 


oU 

°l/2 


£>is(1520) 
511(1535) 
Su(1650) 
£>i S (1675) 
Pis (1680) 
Pis(1720) 


-0.53 1.55 
4.75 1.69 
0.13 1.55 
0.01 2.00 
0.00 1.20 
12.7 1.55 


0.58 1.75 

0.01 2.00 
4.09 1.75 
4.99 1.55 


15.7 1.57 
0.36 1.55 
-0.50 1.55 
0.00 0.00 
0.00 0.00 
0.00 0.00 


13.6 
28.5 
10.1 
0.00 
0.00 
0.00 



Fig. [23J displays our super-global and single-Q 2 fits for the 
D 13 (1520) and F 15 (1680) resonances. The figure demon- 
strates that (I) the helicity non-conserving amplitude A 3 / 2 
dominates for real photons and (II) with increasing val- 
ues of Q 2 , A 3 / 2 drops faster than the helicity conserving 
amplitude A 1 / 2 - As a consequence the asymmetry 



A(Q 2 ) 



Al/2 



A 3/2 I 2 



Al/2 



A m | 2 



(48) 



changes rapidly from values close to —1 to values near +1 
over a small Q 2 range. As is seen in Fig. [24] the asymmetry 
crosses the zero line at Q 2 « 0.5 GeV 2 for the £>i 3 (1520) 
resonance and at Q 2 ps 0.8 GeV 2 for the iq 5 (1680). As a 
comparison, the asymmetry A for the Z\(1232) resonance 
is practically constant over this Q 2 range with a value w 
—0.5. This again shows the special role of the A resonance, 
where the helicity conservation is not observed. 
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6 Conclusion 

Using the world data base of pion photo- and electropro- 
duction and recent data from Bates/MIT, ELSA/Bonn, 
MAMI/Mainz, and Jefferson Lab, we have extracted the 
longitudinal and transverse helicity amplitudes of nucleon 
resonance excitation for all the 4-star resonances below 
W = 2 GeV. For this purpose we have extended our 
unitary isobar model MAID and parameterized the Q 2 
dependence of the transition amplitudes. The comparison 
between such super-global solutions with the correspond- 
ing single-Q 2 fits gives us confidence in the obtained 
helicity couplings for the P 33 (1232), P L1 (1440), 5n(1535), 
£>i3(1520), and Fis(1680) resonances, even though the 
model uncertainty is still quite large, particularly for the 
longitudinal amplitudes. 

For the higher 4-star and all 3-star resonances the 
situation is less clear. This deplorable situation reflects 
the fact that a model-independent analysis requires 
precision data over a large kinematical range. In some 
cases double-polarization experiments will be helpful, 
as has already been shown for pion photoproduction. 
Furthermore, charged pion electroproduction data are 
needed with the same quantity and quality as for neutral 
pions, in order to resolve the ambiguities in the isospin 
structure, in particular for the Sn and S31 resonances. 
While we have mostly discussed the electroproduction 
from proton targets, also the existing neutron data have 
been analyzed. The latter are of course less abundant, 
and moreover no new neutron data have been reported 
over the recent years. Because the isospin symmetry is 
most likely on safe grounds in the resonance region, only 
the electromagnetic neutron couplings with isospin 1/2 
are still lacking. In spite of the discussed problems, we 
have implemented a super-global solution also for the 
neutron amplitudes in MAID07. 

Pion photo- and electroproduction are invaluable tools 
to study the resonance structure of the nucleon. With the 
advent of the new c.w. electron accelerators, new precision 
experiments have afforded a host of new data to unravel 
this structure in the first and second resonance regions. In 
particular, electroproduction has provided new insights in 
the spatial distribution of the nucleon-resonance transi- 
tion densities. However, in order to get the whole picture, 
several challenges remain: 

— Dedicated experiments to investigate the higher en- 
ergy region, which have to include an intense study 
of the polarization degrees of freedom. Experience has 
shown that even the physics of the A (1232) requires 
a full-fledged program to measure the spin observables 
in order to understand the background of the non- 
resonating multipoles. 

— A fresh approach to also determine the excitation spec- 
trum of the neutron. As an example, the comparison 
of the Roper or Pn (1440) helicity amplitudes for pro- 
ton and neutron will shed light on the structure of this 
enigmatic resonance. 



— The open question of the excitation spectrum in the 
third resonance region and above deserves further 
studies in both theory and experiment. This includes 
"missing" and "exotic", e.g., 5-quark resonances as 
well as more mundane second and third resonances in 
a multipole, which show up in a particular analysis and 
not in another one. 

In conclusion we hope that MAID2007, just as other 
approaches based on partial-wave analysis, dynamic 
models, coupled-channels calculations, and dispersion 
theory, will contribute to settle the mentioned issues and 
thus to improve our still somewhat vestigial knowledge of 
the nucleon's resonance structure. 
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Fig. 1. The Eo+ multipole for fp — > 7r°p. Real part: the re- 
sults of MAID98 (dotted line) and MAID2007 without the cusp 
effect (dashed-dotted line) as well as the full MAID2007 cal- 
culation (red solid line). The red dashed line is the imaginary 
part of the full MAID2007 solution. The data points are from 
Refs. rjD](») and [3T](o). 
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Fig. 2. The global solutions of MAID2007 (solid lines) and GWU/SAID [29] (red dashed-dotted lines, solution FA06K) for the 

1/2 3/2 3/2 3/2 

multipoles P E ' + , E ' + , , and as function of the photon lab energy _E 7 in the first resonance region. The blue open 

circles show our single-energy solution. The dashed lines represent our unitarized background contributions in the M\+ and 
E\^ multipoles. 
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Fig. 4. The multipoles n E^ , n M^2 , nE\?, and n M\^ as function of the photon lab energy E-y. Further notation as in Fig. [2] 
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Fig. 6. The global solutions of MAID2007 (solid lines) and GWU/SAID [29] (red dashed-dotted lines, solution FA06K) for the 
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multipoles P E Q ' + , E ' + , M x { , and EA as function of the photon lab energy E 1 in the second and third resonance regions. 



The blue open circles show our single-energy solution. The dashed lines represent our unitarized background given by Eq. ([15 

3/2 

Note that the background for Re M 1+ is out of scale. 
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Fig. 7. The global solutions of MAID2007 (black solid lines) and GWU/SAID 29\ (red dashed-dotted lines) for the multipoles 
p Mj/ 2 , Mj 3 ^ 2 , pEl^ 2 , and V M\^ as function of the photon lab energy Ej in the second and third resonance regions. Further 
notation as in Fig. [6] 
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Fig. 8. The global solutions of MAID2007 (black solid lines) and GWU/SAID [29] (red dashed-dotted lines) for the multipoles 
and as function of the photon lab energy _E 7 in the second and third resonance regions. Further 

notation as in Fig. [6] 
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Fig. 9. The global solutions of MAID2007 (solid lines) and GWU/SAID [29] (red dashed-dotted lines) for the multipoles p e\^ , 

1/2 1 /2 1 II 

pM 2+ , P E 3 _ , and P M 3 _ as function of the photon lab energy _E 7 in the second and third resonance regions. Further notation 
as in Fig. [6] The resonance contribution to the pE^v 2 multipole is very small, and therefore the solid and dashed lines coincide. 
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Fig. 10. The global solutions of MAID2007 (solid lines) and GWU/SAID [29] (red dashed-dotted lines) for the multipoles , 

3/2 3/2 3/2 

M 3 _ , E 3+ , and M 3+ as function of the photon lab energy in the second and third resonance regions. Further notation as 
in Fig. [6] The resonance contribution to the Eg^ 2 multipole is very small, and therefore the solid and dashed lines coincide. 
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Fig. 11. The global solutions of MAID2007 (solid lines) and GWU/SAID [29] (red dashed-dotted lines) for the multipoles 
nE\+ i an d n.A^i+ 2 as function of the photon lab energy E 1 in the second and third resonance regions. Further 
notation as in Fig. [6] 
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Fig. 12. The global solutions of MAID2007 (solid lines) and GWU/SAID [29] (red dashed-dotted lines) for the multipoles 
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notation as in Fig. [6] 



, and nM^ 2 as function of the photon lab energy E 1 in the second and third resonance regions. Further 
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Fig. 13. The global solutions of MAID2007 (solid lines) and GWU/SAID [29] (red dashed-dotted lines) for the multipoles 
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Fig. 14. The contribution of a third Su resonance in the P E Q+ multipole with M_r=1950 MeV, A{=200 MeV, single-pion 
branching ratio /3^=0.4, and helicity amplitude ^1/2=0.028 GeV _1//2 . The solid and dashed lines are our global solutions 
with and without this resonance, respectively. The red dashed-dotted lines represent the global SAID solution. The blue open 
circles and green crosses are the single-energy solutions of MAID2007 and the SAID, respectively. This resonance is included in 
MAID2007. 



28 



D. Drechsel, S. S. Kamalov, L. Tiator: Unitary Isobar Model - MAID2007 




1.0 



+ 


0.5 




E 






i 

o 


0.0 


Si 


CM 


-0.5 






800 900 1000 1100 1200 1300 1400 1500 
E (MeV) 



-1.0 h 
-1.5 

800 900 1000 1100 1200 1300 1400 1500 
E (MeV) 



1/2 - — 

Fig. 15. The contribution of a second P\\ resonance in the V M{'_ multipole with Af fl =1700 MeV, r tot =47 MeV, single-pion 
branching ratio [3^=0.1, and helicity amplitude A\i% = —0.024 GeV _1//2 . The solid and dashed lines are our global solutions 
without and with this resonance, respectively. Further notation as in Fig. 1141 The Pn(1700) is not included in MAID2007. 



D. Drechsel, S. S. Kamalov, L. Tiator: Unitary Isobar Model - MAID2007 



1.2 | i i i i | i i i i | i i i i i i i i i | i i i i | i i i i | i i i i 




12 3 4 5 6 7 

Q 2 (GeV/c) 2 

Fig. 16. The Q 2 dependence of the magnetic form factor 
G* M for the iVZ\(1232) transition divided by 3G D (Q 2 ). The 
solid and dashed blue lines are the results of MAID2007 and 
MAID2003, respectively. The red open triangles represent the 
new JLab data of Ungaro et al. |13j . See Ref. [10] for the no- 
tation of the other data points. 
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Fig. 17. The Q 2 dependence of the ratio Rem at the Zi(1232) resonance. The blue solid and dashed lines are the super-global 
solutions from MAID2007 and MAID2003, respectively. The data points are from Refs. [42] (open square), [43] (solid triangles), 
|47] (open diamond), |48j (cross), [5] (open circle), [62] (asterisks), and [13] (open triangles). The green solid circle at Q 2 = in 
the left panel is from Ref. [4], and the black solid circles in the right panel are obtained by our single-Q 2 analysis of the data 
from Ref. [B]. 


-5 

_ -10 

V 15 f 

\ -20 [ 
" -25 
-30 
-35 




_i i I i i ■ I i ■ ■ I i i i I i i i I i_ 

0.0 0.2 0.4 0.6 0.8 1.0 
\2 



ill! 


1 1 1 1 1 1 1 1 1 1 













» 


; i 




















Fv 




■ . . . i 


J 1 I 1 1 1 I 1 1 L 


. 1 , . . 1 


\ 

A. . • 



1 



Q 2 (GeV/c) 2 



2 3 4 5 6 7 
Q 2 (GeV/c) 2 
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Fig. 19. The Q 2 dependence of the 3 helicity amplitudes for the Zi(1232) resonance, in units 10" 3 GeV~ 1/2 . The solid and 
dashed lines are the MAID2007 and MAID2003 super-global solutions, respectively. The data points are from our sing le-Q 2 fits 
to the 7T° and 7r + CLAS data (red full and blue open circles, see Table [5] for references), Ref. [62] (blue full circles), Ref. [47] 
(black full circles at Q 2 = 0.127 GeV 2 ), and Ref. [2] (green full circles at Q 2 = 0). 
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Fig. 20. The Q 2 dependence of the helicity amplitudes for the Pn(1440) resonance of the proton. Further notation as in Fig. 1191 
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Fig. 21. The Q 2 dependence of the helicity amplitudes for the 5*11 (1535) resonance of the proton. Further notation as in Fig. 1191 
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Fig. 22. The Q 2 dependence of the helicity amplitudes for the Pn(1440) and 5*11(1535) resonances of the proton. The MAID2007 
super-global analysis (solid lines) and the single-Q 2 fits (red full circles with error bars) are compared to the results of Aznau- 
ryan [64] obtained from a similar data set within an isobar model (full triangles) and dispersion theory (open triangles). Further 
notation as in Fig. 1191 
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Fig. 23. The Q 2 dependence of the helicity amplitudes for the Di3(1520) and Fis(1680) resonances of the proton. Notation as 
in Fig. [H 
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Fig. 24. The helicity asymmetry A(Q 2 ) of Eq. (gHJ) for the L>i 3 (1520) and Fi 6 (1680) resonances of the proton. The solid and 
dashed curves are the super-global MAID2007 and MAID2003 solutions, respectively. The data are the results of our single-Q 2 
fits to the CLAS data (red full circles, see Table [9] for references). 



